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1 Introduction 

As was discussed recently [1, 2, 3, 4, 5], chiral perturbation theory [6, 7], the low 
energy effective theory of QCD with Nf light quark flavors, could posses different 
behavior for Nf=2 and Nf=3. As a consequence of vacuum fluctuations of the 
growing number of light quark flavors, the most important order parameters 
of spontaneous chiral symmetry breaking (SBxS), namely the pseudoscalar 
decay constant and the quark condensate in the chiral limit, obey paramagnetic 
inequalities F (N f + 1) < F (Nf) and E(JV> + 1) < E(iV/) [3]. In particular, 
the fluctuations of the sea ss-pairs need not be suppressed due to the relatively 
small value of the s-quark mass m s <AQCD and could bring about a possibly 
significant suppression of Fo(3) and E(3) w.r.t. Nf—2. This should manifest 
itself through the OZI rule violation in the scalar sector as can be seen from 

F Q (2) 2 = F (3) 2 +Wm s B Ll-2fl K + O(m 2 s ) (1) 

E(2) = E(3)(l + ^^Ll - 2fl K - + 0{m% (2) 

where B 2 = E(3)/F (3) 2 , ftp — fJ.p\ mu d ^o- Indeed, L4 and L 6 are the 1/N C 
suppressed LEC's (connected to the scalar mesons), traditionally considered 
negligible. Predictions for L4 and L$ derived from sum rules involving scalars 
[10, 11, 5, 2], calculations on the lattice [12, 13] and NNLO S%PT [14] produce 
numbers significantly different from traditional expectations [7] . Convenient pa- 
rameters relating the order parameters to physical quantities can be introduced: 
Z(N f ) = F (N f ) 2 /F 2 and X(N f ) = 2mZ(N f )/F 2 Ml with m=(m, + m d )/2. 
The large ss vacuum fluctuations could lead to Z(3) < Z(2), A(3) < X(2). 
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Analysis [8] of the K ei decay experimental results for the irir s-wave scattering 
length [9] lead to values X(2)=0.81±0.07, Z(2)=0.89±0.03. ttvt and irK scatter- 
ing data constrain the three flavor parameters much less strictly [1, 4], AT(3)<0.8, 
Z(3)~0.2-0.9, Y=X(3)/Z(3)<1.1, r=m s /m>15. Sum rule approaches [5, 2] 
yield approximately X(2),Z(2)~0.9 and X(3),Z(3)~0.5 at r=25. 

Small X(3) and Z{3) would lead to irregularities of the chiral expansion 
connected to numerical competition of the LO and NLO terms, which could 
be consequently seen as unusually large higher order corrections. An alterna- 
tive approach, dubbed 'Resummed' xPT (R^PT), has been introduced recently 
[1]. It takes this possible scenario into account and is based on the effective 
resummation of the vacuum fluctuation discussed above. The goal of the article 
is to illustrate the 'Resummed' approach on the sector of decay constants and 
to try to use it for theoretical predictions of the 77 decay constant and related 
parameters, including uncertainty estimates. 

The SU(3)l x SU(3)r r\ decay constant in the isospin limit 

ip^Fr, = (0|4h(p)>, (3) 

where are the QCD axial vector currents, can be calculated in SU(3)l x 
SU(3)r xPT without the introduction of the 7/ meson. In the usually investi- 
gated 77-77' mixing sector, the following definitions are used 

iP^Z = (0\A 8 f\ V ,v'). (4) 

As can be seen, the SU(3)l*SU(3)r constant F v is defined identically to F^ 
in the U(3)lxU(3)r framework. A general two angle mixing scheme [15, 16] 

F% = F 8 cosd 8 , F% = F 8 sin7? 8 , F" = -F sin O , F§ = F cost9o (5) 

has been shown to provide better agreement with experimental data and xPT 
predictions [16, 17, 18] than a single mixing angle scenario. Table 1 collects some 
recent two angle phenomenological analyses leading to a value of F^. Older one 
angle mixing scheme results generally provided much lower numbers Fi ~ F n . 

Several recent xPT results can be cited. Standard xPT to 0(p 6 ) [22] gives 
F v /F n = 1 + 0.242 + 0.066 = 1.308. Large N c X PT to NNLO [23] leads to 
F a = lMF n , = -22° and thus = 1.24F W . We build on the 'Resummed' 
XPT result [1] F% = ^(1.651 + 0.036F) (at r=24, remainders neglected). 

2 Decay constants in 'Resummed' %PT 

'Resummed' xPT [1] starts from the same form of the effective Lagrangian as 
the Standard variant (SxPT) [7] . The difference is in the treatment of the chiral 
series, RxPT assumes possible irregularities. Overall convergence to all orders is 
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1.26JV 


-21.2° 


1.17F, 



Table 1: Recent two angle 77-77' analyses leading to a value of 



taken for granted, but only for expansions directly obtained from the generating 
functional. These 'bare' expansions are then dealt with additional caution. 

The first step is to derive a strict chiral expansion fully expressed in terms 
of the original parameters of the effective Lagrangian. In our case we have 



F 2 (l - 4u v - 2 m ) + \&B m{U{r + 2) + L 5 ) + A 



(4) 



(6) 



F'k = F^l-^-3 f x K -^ ri ) + 16B m(U(r + 2) + ^L 5 (r + l)) + A^ K (7) 



^ = *o (1 - + 16B m(L 4 (r + 2) + -L 5 (2r + 1)) + A^ } . 



(8) 



The expansions for the squares of the decay constants are used, as they are di- 
rectly related to two point Green functions obtained from the generating func- 
tional. At this point, the chiral logs fip = m 2 P / 32-k 2 Fq \n{rn 2 P / ix 2 ) contain non- 
physical 0(p 2 ) masses m\ = 2B a rh, m 2 K = B Q rh{\ + r), m 2 = 2/3B rh(l + 2r). 

denote the higher order remainders, not neglected in this approach. 

The second step is the definition of the bare expansion, which usually involves 
changes to the strict form in order to incorporate additional requirements, such 
as physically correct analytical structure. In our case this narrows down to a 
question, whether to replace the original leading order masses inside the chiral 
logarithms with physical ones. In some cases (see [28]) this is a nontrivial 
question, so we will keep both options and evaluate them. 

The next stage is the reparametrization of the unknown LEC's in terms of 
physical observables. In RxPT the leading order ones are left free, only re- 
expressed in terms of more convenient parameters r, Z and X resp. Y . Two 
NLO LEC's are present in our formulae, the equations for F v and F K (6,7) can 
be used for the reparametrization. Note that this is done in a pure algebraic 
way, no additional expansion is made. The final formula for the 77 decay constant 
[1] is then obtained by insertion into (8) 
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The expression is valid to all orders, it's only divided into an explicitly calculated 
part and the unknown higher order remainders A^P . 

The last step consists of the treatment of the remainders. We will use three 
ways to estimate them. The first relies on an assumption about the convergence 
of the chiral series [1, 4] and assumes the typical size of the NNLO remainders 
is |A^| <~ O.lFp. These are added in squares to obtain the final uncertainty. 
The result is a prediction in the sense that a value significantly outside of the 
resulting variance is not compatible with such an assumption about a reasonably 
quick convergence of the chiral expansion. 

Then we try to use information outside core x?T to get a feeling about 
remainder magnitudes. As can be seen, the Rx?T framework is very suitable 
for incorporating such additional sources of information. We collect various 
published estimates for L§ and use them to check the remainder differences 
A^-A^J and A^-A^J. We also use the Generalized %PT Lagrangian [24, 25] 
to get a sense of the magnitude of the higher order corrections 

More details about this procedure will be published elsewhere [28]. 

3 Numerical results 

For the numerical results we use the physical values M 7r = 135MeV, M^=496MeV, 
M^=548MeV, M p =770MeV, F 7r =92.4MeV and F K =113MeV. At first, let us 
investigate the NLO Standard xPT. There are several differences compared to 
the procedure outlined in the previous section. One can use the quadratic form 
of the expansion obtained from the two point Green function or a linearized 
form, as is more usual. For the LEC reparametrization inverted expansions 
for Fq and 2B n rh are used, while r is fixed at r = r 2 = 2Mj c /M 2 — 1 or 
r = f 2 = 3M 2 /2M 2 — 1/2. One then obtains the following formulae 

F n i ,o„ o„ , 8Aff(r-l) . , , , 16jg(r-l) 

where the chiral logs contain physical masses only fip = Mp/32ir 2 F% \n(Mp / /j, 2 ) . 

As for LI, we opted to use the published values L r b {M p ) = (1.4 ± 0.5).10 -3 
(C(p 4 ) fit [7, 26]) and L$(M P ) = (0.65 ± 0.12). lO" 3 (0(p 6 ) fit [27]). 

All these possibilities differ merely in redefinitions of the usually neglected 
remainders. Table 2 shows that it might be worth to spend the additional effort 
to bring the higher order uncertainties explicitly under control. Numerically, the 
sensitivity to the change in L r 5 is in the range AF^/F^ = (0.11 — 0.14) Ai 5 .10 3 . 
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0{p A ) L s , r=r 2 


0(p 6 )L 5 , r=r 2 


0(p 4 ) L 5 , r=r 2 


0(p 6 )L 5 , r=f 2 


Fr, 


1.31±0.07 


1.21±0.02 


1.29±0.07 


1.19±0.02 


VF* 


1.27±0.06 


1.19±0.01 


1.25±0.05 


1.17±0.01 



Table 2: Various NLO S%PT results for the rj decay constant in F v units. 



Proceeding to RxPT, we generally investigate a standard and a low r sce- 
nario r ^15-25 and vary Y in the range 0-1.6. Keep in mind, though, that the 
7T7T and irK scattering analyses [1, 4] suggest Y < 1.1. 

Let us first neglect the remainders and have a look on the dependence of 
the explicitly calculated part on the free parameters Y, r and the treatment of 
chiral logs. As can be seen from Fig.l, the dependence on both is very small. 
For physical masses inside the logs one gets F%/F% = 1.661 — 0.011F+0.002Fr. 
The decay constant sector might thus be insensitive to the particular scenario of 
SBxS and more information is needed to extract the values of the parameters. 
Neglecting these weak dependencies one gets the sensitivity on the remain- 



ders as AF V /F n = 1.5 • 10-y((3A^ 4) ) 2 + (4A^) 2 + (A^) 2 ) 




1.26 



Figure 1: F n in RxPT, remainders neglected. 
Chiral logs: solid - physical masses, dashed - 
0(p 2 ) masses. Dark: r = 25, light: r = 15 



Applying the 10% uncertainty 
remainder estimate, the fol- 
lowing all order approxima- 
tion is obtained 

F v = (1.3±0.1)F W . (12) 

All phcnomcnological and 
theoretical results cited in the 
introduction fall in or very 
close to this range and are 
thus compatible with a reason- 
able convergence of chiral se- 
ries. However, the mentioned 
one mixing angle scheme re- 
sults are significantly outside. 



The difference F\-F 'I 



mate on the remainder difference A^P- 



(6,7) depends only on L\. This yields an order esti- 



A^' if independent information on L\ 



can be gathered. Several estimates for L r b beyond 0(p 4 ) xPT are available: 

- SxPT 0(p 6 ) fit: L r 5 (M p ) - (0.5 - 1.0) • 10~ 3 [27, 22] 

- Resonance saturation: LI ~ (1.6 — 2.1) • 10~ 3 [11] 

- QCD sum rules: L r 5 (M p ) > 1.0 • 10~ 3 [2, 5] 

- xPT on lattice: L r 5 - 1.8 - 2.2 • 10~ 3 [12, 13] 

The result of varying L£(M p ) in the range (0.5 — 2). 10 -3 can be seen in Fig. 2. 
0(p 2 ) masses were kept inside logarithms, physical ones make the remainder 
estimate somewhat larger. The estimate is compatible with small remainders. 
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Figure 2: A^-A^j estimate for L r 5 (M p ) ~ (0.5 - 2).1(T 3 (dark band). 
Upper bound corresponds to low values of L\. Left: r=15, right: r=25. 
Light: expected uncertainty ±0.12F|- from the 10% uncertainty estimate. 
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Figure 3: A^j-A^j estimate for L r 5 (M p ) ~ (0.5 - 2).10- 3 (dark band). 
Upper bound corresponds to low values of L\. Left: r=15, right: r=25. 
Light: expected uncertainty ±0.11F 2 from the 10% uncertainty estimate. 



We can also utilize the information about the difference F^-F% (6,8). If we 
use the latest phenomenological result F* ~ 1.38F,r [18] as an input, an estimate 

of Ap^-Ap^ is obtained. It should be stressed that older results produced 
lower values, so this should be taken as a preliminary look on the possible 
consequences if such a higher value of F v was confirmed. We don't make a full 
statistical analysis, only provide some first feelings where it could lead to. Keep 
in mind |A^| ~ 0.1 Fp and Y < 1.1 as suggestions following from [1, 4]. These 
assumptions hint the following consequences, demonstrated in Fig. 3 

- r-15 and A^-A^O.2^ 2 implies Y>1 or L r 5 (M p )>2. 10~ 3 . 

- r-25 and A^-A^O^F 2 implies F>0.5 or L£(M p )>2.10- 3 . 

- L r 5 (M p ) <1.10" 3 and r-25 implies Y>1.2, A^-A^O.2^ 2 . 

The remainder estimate using the Generalized xPT Lagrangian [28] provides 

3A^ ) -4A^]+A^ ) = 2^F 2 m 2 (r-l) 2 +8^ J F l2 TO 2 (r 2 + l)-4B^( M ) J F l2 m 2 (r-l) 2 
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Figure 4: GxPT remainder estimate for r—25. Dark: Z—0.9, light: Z=0.5. 
Left: simple variation of scale, solid: [i = lGcV, dashed: [i = M p . 
Right: the LEC estimate described in the text, solid: error bars \i = lGeV/M p , 
dashed: central values /x = M p . 



We use two ways to estimate the unknown G^PT LEC's. The first is the 
usual simple variation of scale, the constants are set to zero at two different scales 
and the sensitivity is checked. The second assumes a probabilistic distribution 
of possible values depending on scale variation 

= "W 1 " 1 ^ 1 ' ™ = 0± fef h O- (14 > 

These are then added in squares. Note that the insensitivity in the first case 
assures independence on where the central value is chosen in the latter one. 

The results for r=25 can be seen in Fig. 4, low values of r do not change 
the overall picture. However, of the four G%PT LEC's present in our case only 
two depend on scale, which is hardly a good statistical ensemble. There is no 
indication of large higher order corrections nevertheless. 

4 Summary 

We have studied the case of pseudoscalar decay constants in the 'Resummed' 
%PT framework and tried to obtain an estimate for the rj decay constant and 
related parameters. We used several ways to get a feeling about the effect of 
higher order remainders. 
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